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We apply in this paper the statefinder parameters to the interacting phantom energy with dark 
matter. There are two kinds of scaling solutions in this model. It is found that the evolving 
trajectories of these two scaling solutions in the statefinder parameter plane are quite different, and 
that are also different from the statefinder diagnostic of other dark energy models. 

PACS numbers: 98.80.-k, 98.80.Es 

Cosmic observations indicate that our universe is undergoing an accelerated expansion and the dominated com- 
ponent of the present universe is dark energy 

GL H SHI!- The Wilkinson Microwave Anisotropy Probe (WMAP) 
satellite experiment tells us that dark energy dark matter, and the usual baryonic matter occupy about 73%, 23%, and 
4% of the total energy of the universe, respectively. The accelerated expansion of the present universe is attributed to 
the dark energy whose essence is quite unusual and there is no justification for assuming that it resembles known forms 
of matter or energy. Candidates for dark energy have been widely studied and focus on the cosmological constant A 
[1,0] with W = — 1, a dynamically evolving scalar field (quintessence) [1, Q with W > — 1 and phantom with 
W < — 1. Recently, a study of high-Z (Z is redshift) SNe la [ll[ find that the equation of state of dark energy has a 
99% probability of being VV < — 1 if no priors are placed on Q^, . When these SNe results are combined with CMB 
and 2dFGRS the 95% confidence limits on an unevolving equation of state are —1.46 < W < —0.78 d, HH which is 
consistent with estimates made by other groups [4, 5]. In order to obtain W < — 1, phantom field with reverse sign 
in its dynamical term may be a simplest way and can be regarded as one of interesting possibilities describing dark 
energy However, the other physical properties of phantom energy are rather weird, as they include violation 

of the dominant-energy condition, naive superluminal sound speed and increasing energy density with time. The 
latter property ultimately leads to unwanted future singularity called big rip which had been considered in [l3| . This 
singularity is characterized by the divergence of the scale factor in a finite time in future [3] • 

Many authors have discussed various kinds of phantom field models to avoid the cosmic doomsday (Til ] , which require 
a special class of phantom field potentials with a local maximum. Moreover, the energy density of the phantom field 
increases with time, while the energy density of the matter fluid decreases as the universe expands. Why are the 
energy density of dark matter and the phantom energy density of the same order just at the present epoch? This 
coincidence problem becomes more difficult to solve in the phantom model without the suitable interaction 16]. But 
Guo et al. in Ref. [ij], [l|| proposed a suitable interaction in the phantom field model, and the coincidence problem 
can be avoided. Moreover in Ref. [HI], the universe also avoids the big rip. In Ref. [HI, considering a universe model 
which contains phantom field 4> and the dark matter pdm- The Friedmann equation in a spatially flat FRW metric 
can be written as follows: 



H 2 = —{p^+ Pdm), (1) 



and 



H=~(p dm -4> 2 ), (2) 

where k 2 = 8irG, pdm is the energy density of the dark matter, and the dark matter possesses the equation of state 
Pdm = 0. The energy density and pressure of the phantom field <f> are p^ and P^, respectively, 

P<b = A^ + V{4>), (3) 

p = -^ 2 -n# (4) 
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where V(<p) is the phantom field potential, V((j)) — Vq exp(— Ak</>). Here, we postulate that the two component p$ and 
Pdmi interact through the interaction term Q according to 

p <j> + 3H(p < f > +P <t> ) = -Q, (5) 



Pdm + ZH(p dm + P dm ) = Q. (6) 

Here, the interaction term has the specific formation 

Q = 3cH(p ( / > + pdm)- (7) 
where c is a dimensionless parameter denoting the transfer strength. We define the following dimensionless variables 



KVV KjPdm ,_, 

y=^77' z = Krr ■ ( 8 ) 



y/GH' V3H' ' V3H 

Thus the fractional densities of p^ and pd m are 

n <t> = -x 2 + y 2 , n dm = z 2 . (9) 

The evolution equations (J3j) and (6|) can be written as the following set of equations: 

x = -3z(l + x 2 - \z 2 \cx~ 2 ) - ^Ay 2 , (10) 

V = -3 y (x 2 + ^Ax-iz 2 ), (11) 

z = -3z(i + x 2 -iz 2 -icz- 2 ), (12) 

where the prime denotes a derivative with respect to the logarithm of the scale factor, N = In a. From the definitions 
of these new variables, we find that the equation of state of phantom is = (— x 2 — y 2 )/(—x 2 + y 2 ). So we get 
the effective parameter of equation-of-state of the phantom is 

n ff - " x2 ;fr - d3) 

v —x z + y z 

and the effective equation of state for the total cosmic fluid is 

W e// = -x 2 - y 2 . (14) 

Guo et al. in Ref. [18| had concluded that in the case of the interaction ([7]), there exist two stable scaling solutions, 
the climbing-up scaling solution and the rolling-down scaling solution. In this model the universe evolves from a 
matter-dominated phase to a scaling solution, which is characterized by a constant ratio of the energy densities of 
the dark matter and the phantom field, which may give a phenomenological solution of the coincidence problem. In 
the climbing-up case, the phantom field initially climbs up, the effective equation of state W e // realizes a transition 
from Weff > — 1 to W e // < — 1, and the universe ends with a big rip. In the rolling-down case, the phantom field 
initially rolls down, the effective equation of state tends to above —1 and realizes a transition from W e // < — 1 to 
W e // > —1, in this case the cosmic doomsday is avoided and the universe accelerates forever. There have been so 
many models that proposed to explain the cosmic acceleration and solve, or at least alleviate, the coincidence problem 
(9l UM. [20l| . In order to discriminate this interacting phantom model from others, and differentiate these two scaling 
solutions further, we refer to a cosmological diagnostic pair {r, s} called statefinder which is introduced by Sahni et 
al. in [2l| and defined as 

a r — 1 , , 

(15) 



~ aH 3 ' ~ 3(^-1/2)' 

Here q is the deceleration parameter. The statefinder is a "geometrical" diagnostic in the sense that it depends on the 
expansion factor and hence on the metric describing space-time. Since different cosmological models involving dark 
energy exhibit qualitatively different evolution trajectories in the s—r plane, this statefinder diagnostic can differentiate 
various kinds of dark energy models. For the spatially flat LCDM cosmological model, the statefinder parameters 
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FIG. 1: The left figure is s — r diagram of the climbing-up scaling solution with xo < 0. The right figure is q — r diagram of the 
climbing-up scaling solution with xo < 0. The curves evolve in the variable interval N £ [—2,20]. Selected curves for A = 1, 
c = 0, 0.1, 0.2 and 0.3 respectively. Dots locate the current values of the statefinder parameters. 



correspond to a fixed point {r = 1, s = 0}. S far, some models including the cosmological constant, quintessence, 
phanto m, q uintom, the Chaplygin gas, braneworld models, holographic models, interacting and coupling dark energy 
models dSHH, HI], have been successfully differentiated. For example, the quintessence model with inverse power 
law potential, the phantom model with power law potential and the Chaplygin gas model all tend to approach the 
LCDM fixed point, but for quintessence and phantom models the trajectories lie in the regions s > 0,r < 1 while for 
Chaplygin gas model the trajectories lie in the regions s < 0, r > 1. In this paper we apply the statefinder diagnostic 
to the coupled phantom model. To begin with, we use another form of statefinder parameters in terms of the total 
energy density p and the total pressure p in the universe: 

r = 1+ 9(p + p)g s= iP±At. (16) 
2p p VP 

Since the total energy of the universe is conserved, we have p — —3H(p + p). Then making use of p^ = —3H(1 + 

r = 1 - ' d -W'^ + ^(1 + W e / f )a$, (17) 

where = -j^- and the deceleration parameter is also given 

g = i(l + 3W*n ) (19) 

In the following we will discuss the statefinder for two scaling solutions with different conditions in Ref 18], and 
investigate how the interaction between phantom and dark matter influences the evolution of the universe. Firstly, 
we discuss the scaling solution B in the case of the interaction form ([7]), which exists for ig < (xb is the initially 
velocity of the phantom field) and < c < /( ~ A ~/^ + -) {fix) = x(2x + V6A/3)(1 - %/6xA/3) is defined as a cubic 
function) and corresponds to a climbing-up phantom field. In this scenario, the universe ends with a big rip. In Fig. 
[T]we show the time evolution of the statefinder pairs {r, s} and {r, q} for the climbing-up scaling solution. The plot is 
for variable interval N e [—2, 20], and the selected evolution trajectories of r(s) and r(q) correspond to A = 1, c = 0, 
c = 0.1, 0.2 and 0.3 respectively, c = represents no interaction between phantom field and dark matter, c = 0.1, 0.2 
and 0.3 represent the different transfer strength of interaction between DE and DM. We see clearly that the curve will 
pass through the LCDM fixed point when there is no coupling (c = 0) between phantom field and dark matter, while 
the distant from the curves to LCDM scenario is considerable when exists the interaction between phantom field and 
dark matter. The statefinder pair {r, s} lies in the regions s < 0, r > 1, which is different from other quintessence and 
phantom model. In Fig. [5] we plot the evolution trajectories of the statefinder parameters versus redshift diagram 
of the climbing-up scaling solution with xo < 0. In the plots of the — ln(l + z) — r and — ln(l + z) — s, we can 
see clearly that the interaction between phantom field and dark matter causes the big deviation between statefinder 
parameters and the LCDM scenario. 

Next we discuss another scaling solution C in the case of the interaction form (J7J) in [l8| , which exists for xq > 
(x c is also the initially velocity of the phantom field) and < c < min{f( ~ A+ ^ 2+6 ), /( ~ A+ 2 ^ +12 )} , which 
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FIG. 2: The figures are the statefinder parameters versus redshift diagram of the climbing-up scaling solution with xo < 0. 
The left figure is — ln(l + z) — r diagram in the variable interval N £ [—2, 20]. The middle figure is — ln(f + z) — s diagram 
in the variable interval N € [—2, 20]. The right figure is — ln(l + z) — q diagram in the variable interval N G [—2, 4]. Selected 
curves for X = 1, c = 0, 0.1, 0.2 and 0.3 respectively. Dots locate the current values of the statefinder parameters. 

corresponds to a rolling-down phantom field where universe avoids the big rip. In Fig. [3] we show the time evolution 
of the statefinder pairs {r, s} and {r, q} for the rolling-down scaling solution. The plot is also for variable interval 
N € [—2,20], and the corresponding parameters are also A = 1, c = 0, c = 0.1, 0.2 and 0.3 respectively. We can see 
that the trajectories of this case will pass through LCDM fixed point. The statefinder pair {r, s} lie in the region s > 0, 
r < 1 which is different from that in Fig. [TJ In Fig. Q]we plot the evolution trajectories of the statefinder parameters 
versus redshift diagram of the rolling-down scaling solution with xq > 0. In the diagrams of the — ln(l + z) — r and 
—ln(l + z) — s, we can see clearly that the interaction between phantom field and dark matter causes statefinder 
parameters to approach the LCDM scenario in the future. In Fig. [FJ the evolution of the effective equation of state 
of the total cosmic fluid with different initial velocity have been shown. In the left figure, the phantom field initially 
climbs up, and the effective equation of state W e // tends to below —1 and realizes a transition from W e // > — 1 
Weff < — 1, thus the universe ends with big rip. In the right figure, the phantom field initially rolls down, the effective 
equation of state W e // tends to above —1 and realizes a transition from W e // < — 1 to W e // > — 1- In this case the 
universe avoids the big rip. 

The scaling solutions B and C in [l8[ with different initially velocity of the phantom field cause the different 
evolution of the universe. We apply a statefinder analysis to this model, and contrast the scenario with coupling 
between phantom field and dark matter and the scenario without interaction between these two components. The 
difference can be found in Fig. [T]-Fig. [5] (i) The region of the statefinder pair {r, s} is different between Fig. [1] 
and Fig. [3J in Fig. [3] s > and r < 1 , while in Fig. [1] s < and r > 1, which differs from the quintessence and 
phantom models, (ii) The influence of the interaction form is different in the evolution of the universe, which can be 
seen in the statefinder trajectories. In Fig. [2] we can see that the interaction form causes the big deviation between 
statefinder parameters and the LCDM scenario. In Fig. 3] the interaction causes statefinder parameters to approach 
the LCDM scenario in the future. In Fig. the effective equation of state W e // realizes a transition from W e // > —1 
to W e // < — 1 and the universe ends with big rip if the phantom field initially climbs up. Contrarily, the effective 
equation of state transits from W e // < — 1 to W e // > — 1 and the universe avoids the big rip if the phantom field 
initially rolls down. So, through the statefinder diagnostic, we not only see the influence of the interaction form on 
the evolution of universe with different initial velocity, but also contrast the difference between the scenario with and 
without the interaction. 

In summary, we study the statefinder parameters of the coupled phantom model. We analyze two cases of scaling 
solutions — the climbing-up scaling solution and rolling-down scaling solution, and contrast the difference between the 
scenario with and without the interaction. It is found that the evolving trajectories of these two scaling solutions in 
the s — r and q — r plane is quite different, and which is also different from the statefinder diagnostic of other dark 
energy models. We hope that the future high precision observation will be capable of determining these statefinder 
parameters and consequently shed light on the nature of dark energy. 
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FIG. 3: The left figure is s — r diagram of the rolling-down scaling solution with xq > 0. The right figure is q — r diagram 
of the rolling-down scaling solution with xo > 0. The curves evolve in the variable interval TV 6 [—2,20]. Selected curves for 
A = 1, c = 0, 0.1, 0.2 and 0.3 respectively. Dots locate the current values of the statefinder parameters. 
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FIG. 4: The figures are the statefinder parameters versus redshift diagram of the rolling-down scaling solution with xo > 0. 
The left figure is — ln(l + z) — r diagram in the variable interval TV £ [—2, 20]. The middle figure is — ln(l + z) — a diagram 
in the variable interval TV £ [—2, 20]. The right figure is — ln(l + z) — q diagram in the variable interval TV G [—2, 4]. Selected 
curves for A = 1, c = 0, 0.1, 0.2 and 0.3 respectively. Dots locate the current values of the statefinder parameters. 




FIG. 5: The figures are the effective equation of state for the total cosmic fluid w e ff versus redshift. The left figure is — ln(l + z) 
— w e f f diagram of the climbing-up scaling solution with xq < 0. The right figure is — ln(l + z) — w e f f diagram of the rolling- 
down scaling solution with xo > 0. The curves evolve in the variable interval TV 6 [—2, 4]. Selected curves for A = 1, c = 0, 0.1 
and 0.2 respectively. 
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